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Abstract 

We analyse a generalisation of general relativity that incorporates 
a cosmic time-variation of the velocity of light in vacuum, c, and 
the Newtonian gravitation 'constant', G, proposed by Albrecht and 
Maguejo. We find exact solutions for Priedmann universes and de- 
termine the rate of variation of c required to solve the flatness and 
classical cosmological constant problems. Potential problems with 
this approach to the resolution of the flatness and classical cosmolog- 
ical constant problems are highlighted. Reformulations are suggested 
which give the theory a more desirable limit as a theory of varying G 
in the limit of constant c and its relationship to theories with varying 
electron charge and constant c are discussed. 

PACS Numbers: 98.80.Cq, 98.80.-k, 95.30.Sf 



1 Introduction 

Albrecht and Magueijo flj have recently drawn attention to the possible cos- 
mological consequences of time variations in c, the velocity of light in vacuum. 
In particular, there are mathematical formulations of this possibility which 
offer new ways of resolving the horizon and flatness problems of cosmology, 
distinct from their resolutions in the context of the standard inflationary 
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universe theory [[J or the pre-Big-Bang scenario of low-energy string theory 
||. Moreover, in contrast to the case of the inflationary universe, varying 
c may provide an explanation for the relative smallness of the cosmological 
constant today. The study presented in ref . |I[ envisages a sudden fall in the 
speed of light, precipitated by a phase transition or some shift in the values 
of the fundamental constants, and explores the general consequences that 
might follow from a sufficiently large change. 

In this paper we explore the formulation of ref. |l] in further detail, show- 
ing that one can solve the cosmological field equations that define it for gen- 
eral power-law variations of c and G. This allows us to determine the rate 
and sense of the changes required in c if the flatness, horizon, and cosmolog- 
ical constant problems are to be solved. We find that these questions have 
simple answers that generalise the conditions for an inflationary resolution 
of the flatness and horizon problems in standard generalised inflation. We 
will also explore some of the consequences of the varying c model for other 
cosmological problems, like isotropy, and its implications for the evolution of 
black hole horizon. We will then look more critically at the formulation of a 
varying-c theory of gravitation that is proposed in ref. [p] , and suggest some 
ways in which it might be improved. In particular, since the proposed formu- 
lation also includes the possibility of varying gravitation 'constant', G(t), we 
consider the relation between this part of the theory and other well-defined 
scalar-tensor gravity theories which incorporate varying G. 

There have been several high-precision tests of any possible space and 
time variations of the fine structure constant at different times in the past, 
@], f|. Although any theory which admits c variation allows us to com- 
pare its predictions with the observational constraints on time variations 
of the fine structure constant, these observational limits arise from relatively 
low-redshift observations or from primordial nucleosynthesis and are not nec- 
essarily of relevance to the varying-c theory under examination here. The 
period of c-variation is expected to be confined to the very early universe 
and its observational consequences are likely to be manifested most sensi- 
tively through the spectrum of any microwave background fluctuations that 
are created. 

2 The Albrecht-Magueijo Model 

If we begin with the field equations of general relativity 
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Gab — g^A — ^ 4 T ab (1) 

then we cannot simply 'write in' variations of constants like G or c as we 
could in a non-geometrical theory like Newtonian gravity ||. For example, 
if we wished to allow G, say, to vary in time then we encounter a consis- 
tency problem if the energy-momentum tensor T a \, and its vanishing covari- 
ant divergence retains its usual form and physical meaning as conservation 
of energy and momentum. The covariant divergence of the left-hand side of 
([]]) vanishes and if the conservation of energy and momentum is expressed 
by the vanishing divergence of T ao in the usual way then G must be con- 
stant. Thus scalar-tensor gravity theories of the Jordan-Brans-Dicke (JBD) 
sort must derive the variation of G from that of a dynamical scalar field 
which contributes it own energy-momentum tensor to the right-hand side of 
the field equations. An alternative way of incorporating a varying G or c 
with a minimum of change to the underlying theory of gravity is to allow the 
conventional energy-momentum tensor to have a non-vanishing divergence. 
Thus the usual energy and momentum will not be conserved. Equivalently, 
one might interpret this as simply changing the definition of the conserved 
quantity. Thus, since the divergence left-hand side of (|l|) must vanish we 
require that the divergence of G(x)c(x) _4 T afe vanishes. 

When discussing any theory in which there is a varying of some tradi- 
tional 'constant' of Nature, it is important to recognise that invariant opera- 
tional meaning can only be attached to space-time variations of dimensionless 
constants. Variations in dimensional constants can always be transformed 
away by a suitable choice of coordinate frame (change of units); this is dis- 
cussed in detail in ref. [[TTJ. We will be interested in theories in which a 
time- variation of the fine structure constant is represented through a vari- 
ation in c. An example (not unique) of a representation as a theory with 
varying electron charge, e, but constant c has been given by Bekenstein, ||, 
and another by Barrow and Magueijo, [|13|]. In varying speed of light (VSL) 
theories a varying a is interpreted as c oc h oc or l l 2 and e is constant, Lorentz 
invariance is broken, and so by construction there is a preferred frame for 
the formulation of the physical laws. In the minimally coupled theory one 
then simply replaces c by a field in this preferred frame. Hence, the action 
remains as 
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with ip(x^) = c 4 . The dynamical variables are the metric g^, any matter 
field variables contained in the matter lagrangian Cm, and the scalar field ip 
itself. The Riemann tensor (and the Ricci scalar) is to be computed from 
g^ v at constant ip in the usual way. This can only be true in one frame: 
additional terms in d^ip must be present in other frames; see refs. fpTII — fp~3|] for 
more detailed discussion. 

Varying the action with respect to the metric and ignoring surface terms 
leads to 

G/_ w — g^i/A = Tpv (3) 

Therefore, Einstein's equations do not acquire new terms in the preferred 
frame. Minimal coupling at the level of Einstein's equations is at the heart 
of the model's ability to solve the cosmological problems. It requires of 
any action-principle formulation that the contribution C$ must not contain 
the metric explicitly, and so does not contribute to the energy-momentum 
tensor. 

Albrecht and Magueijo (AM) propose that a time-variable c should not 
introduce changes in the curvature terms in Einstein's equations in the cosmo- 
logical frame and that Einstein's equations must still hold. Thus, c changes in 
the local Lorentzian frames associated with the cosmological expansion and is 
a special-relativistic effect. The resulting theory is clearly not covariant and 
so implementation of this idea requires a specific choice of time coordinate. 
Choosing that specific time to be comoving proper time, and assuming the 
universe is spatially homogeneous and isotropic, so that there are no spatial 
variations in c or G, leads to the requirement that the Friedmann equations 
still retain their form with c(t) and G(t) varying. Thus the expansion scale 
factor obeys the equations 



a 2 8nG(t)p Kc 2 {t) 



a 2 3 a 2 



(4) 



4irG(t), 3p . 

^iP+^) a ( 5 ) 



c 2 (t) 

where p and p are the density and pressure of the matter, respectively, and 
K is the metric curvature parameter. By differentiating (Q) with respect to 
time and substituting in (||), we find the generalised conservation equation 
incorporating possible time variations in c(t) and G(t), 
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The cosmological constant is a further feature of interest in these theories. 
In general relativity we know that there are two equivalent ways of introduc- 
ing a cosmological constant into the field equations and their cosmological 
solutions. An explicit constant A-term can be added to the Einstein- Hilbert 
gravitational lagrangian. Alternatively, a matter field can be introduced 
which describes a perfect fluid with a stress satisfying p = —pc 2 . In theories 
other than general relativity these two prescriptions need not lead to the 
same contribution to the field equations and their solutions. For example, in 
Brans-Dicke theory or other scalar-tensor theories one can see explicitly that 
they contribute different terms to the cosmological evolution equations (see 
Barrow and Maeda ) . The matter lagrangian term Cm will in general also 
contain large quantum contributions, so Cm = £m [class] + Cm [quantum]. 

If we examine the action given in eq. (|2|) then we see that there are two 
ways in which a cosmological constant can arise in a VSL theory. In our 
discussion we shall model the explicit A term of geometrical origin in (0), 
or (0). This can still be interpreted as a matter field obeying the equation 
of state pa = —p^c 2 , as would be contributed by any slowly changing self- 
interacting scalar field. We do not address the quantum gravitational or 



quantum field theoretic versions of the problem discussed in refs. ||T2| . |8[ 
or solve the problem associated with the presence of quantum contributions 
described by any constant term in Cu\quantum\. 

If we wish to incorporate such a cosmological constant term, A, (which 
we shall assume to be constant) then we can define a vacuum stress obeying 
an equation of state 

Pa = -PaC 2 , (7) 

where 

Ac 2 

Then, replacing p by p + p A in (||) , we have the generalisation 

a, Ps . G 3Kcc 

p + 3-(p +?)+ „ A = -„- + _ (9) 



5 



We shall assume that the remaining matter obeys an equation of state of 
the form 

p = (7 - l)pc 2 (t) (10) 

where 7 is a constant. 

3 The Flatness and Lambda Problems 

If we specialise the standard situation in general relativity in which we set 
G = = c then the conservation equation, ([]), gives p oc a -37 and the 
8irGp/3 term will dominate the curvature term, Kc 2 a~ 2 at large a so long as 
the matter stress obeys p + 3p/c 2 < 0, p + p/c 2 > 0, that is, if 

0<7<jj (11) 

This is what we shall mean by the flatness problem. Since the scale factor 
then evolves as a(t) oc t 2//37 if 7 > or exp{H t}, H constant if 7 = 0, we 
see that it also grows faster than the proper size of the particle horizon scale 
(oc t) so long as 7 < 2/3. Thus, a sufficiently long period of evolution in the 
early universe during which the expansion is dominated by a gravitationally 
repulsive stress with p + 3p/c 2 < can solve the flatness and horizon prob- 
lems. Such a period of accelerated expansion is called 'inflation'. However, 
we notice that if a constant vacuum stress, p A , associated with a non-zero 
cosmological constant is added to the right-hand side of (H) then, in order to 
explain why it does not totally dominate the 87rGp/3 term at large a(t), we 
would need a period during which the early universe was dominated by an 
extreme fluid with 7 < 0, that is p + p/c 2 < 0. This is problematic because 
it leads to contraction a < and apparent instabilities of the vacuum and 
of flat space-time. This is what we shall mean by the cosmological constant 
(or lambda) problem. Other quantum cosmological attempts to solve this 
cosmological constant problem have so far proved unsuccessful ||. 

We see that the inflationary solution of the flatness problem relies on a 
change in the evolutionary behaviour of the matter term in the Friedmann 
equation which allows it to dominate the curvature term at large a(t) and 
transform the zero-curvature solution into the late-time attractor. In con- 
trast, the varying-c model of AM will provide a solution of the flatness and 
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lambda problems by introducing a variation that increases the rate of fall- 
off of the curvature (Kc 2 a~ 2 ) and cosmological constant (p A ) terms in the 
Friedmann equation with respect to the 8nGp/3 term at large a(t). This is 
possible for all values of 7 and does not necessarily require the existence of 
a period during which p + 3p/c 2 is negative. 



4 Exact solutions with Varying G(t) and c(t) 

We need to obtain solutions of eqns. (||) and (|9|) in order to evaluate the 
effects of varying G and c on the expansion dynamics. Let us first consider 
the implications for the flatness and horizon problems alone by setting A = 0. 



4.1 The Flatness Problem 

In order to solve eq. (|j) we assume that the rate of variation of c is propor- 
tional to the expansion rate of the universe; that is, 

c {t) — Coa n ; c , n constants. (12) 

We shall see that we do not need to prescribe a form for the time variation 
ofG(t). 

From eqs. (|S|) and (|TUD, we have 

(pa^y G 3Kcc 



pa?~i G AirGpa 2 ' 

hence 

(Gpa 3 " 1 )- 3Kcc 
Gpa^ ~ AixGpa 2 ' 

and so, using (|T2|), we obtain 



^. = 6Kc^na a 

Aix 

Integrating, we have the exact solutions 

3Kr 2 r>n 2n+3 ^<~ 2 

G(t)pa 3 ^ = + B, if 2n + 3 7 + 2, (14) 

V W 4vr(2n + 37-2) ' ' T ' K ' 
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G(t)pa 3 ^ = °- + B, if 2n + 3 7 = 2, (15) 

47T 

with a constant of integration, B > 0. This holds for arbitrary varia- 
tions. Analogous solutions can be found from (|T3|) for other possible varia- 
tions of c with respect to a(t) or t but we shall confine our attention to the 



simple case of eq. (|l"2"D. 



Returning to the Friedmann eqn., (|J), we see that the condition for the 
solution of the flatness problem is the same as that which would hold if G 
were constant. The only impact of G(t) is to multiply p, so we have in (|]), 
(writing only the 2n + 37 7^ 2 case explicitly) 



4 = l K ^ n f n * +B'a- 3 i -Kc 2 Q a 2n - 2 
a 2 (2n + 37-2) 

- = ^a- 3 ^+ ° 1 5' constant, 

a 2 (2n + 3 7 — 2) 



(16) 



and the B' term dominates the curvature (K) terms at large a, so the flatness 
problem can be solved so long as 

™<^(2-3 7 ), (17) 



where we have noted from ( 15)) that the flatness problem is also resolved in 



the 2n + 3 7 = 2 case. We see that this condition is a straightforward general- 
isation of the inflationary condition (|ll]) that obtains when c is constant. We 
confirm that the variation in G does not enter this condition in a significant 
way. Also, the sign of 37 — 2 determines the overall sign of the curvature 
term in the Friedmann equation. However, unlike in the constant c case, it 
is possible to resolve the flatness problem without appeal to a matter source 
with 3 7 — 2 < 0. In particular, we see that in a radiation-dominated universe 
(7 = 4/3), the flatness problem is soluble if c(t) falls sufficiently rapidly as 
the universe expands, with 

n<-l, (18) 
while for dust (7 = 1) we require a slower decrease, with 

n<-\. (19) 
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4.2 The Horizon Problem 

The variation of the expansion scale factor at large a(t) approaches 

a(t) oc t 2/37 (20) 

when the curvature term becomes negligible, as in the situation when c is 
constant. The proper distance to the horizon increases as 

4 oc ct oc a n t oc t (Wn)/3 ^ (21) 

and so we see that the scale factor can grow as fast, or faster, than dh as t 
increases if 

2 > 2n + 37 

This is the condition to solve the horizon problem and is identical to that 
for the solution of the flatness problem, just as in the situation with constant 
c. 



4.3 The Lambda Problem 

If we consider the case with a non-zero A term, that is contributed by a stress 
obeying p\ = —p A c 2 , then we need to solve (g) with c(t) varying as (|I2"D. In 
this case we need to assume a form for the variation of G with a or t, and we 
shall again assume that its rate of variation is proportional to the expansion 
rate of the universe, with 

G(t) = G a q (22) 

where Go and q are constants. We assume that A is constant. From eq. (^) 
we have 

Ac 2 (2n - g)a 2 "- f '- 1 d 
Pa ~ 8^ 

Now we look for solutions of @ with 

pa 37 = Da- f + Ea- g (23) 

where B, D,and E are constants. 
This requires 
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/ = g-2n-3 7 (24) 

g = 2 + g-2n-3 7 = 2 + / (25) 

D = Acl(q ~ 2U) (26) 

8vrG (2n + 3 7 ) 1 ' 

3Kc 2 n 

47rG (2n + 3 7 - 2) y } 

where we only write down the solution for the cases where 2n + 3 7 ^ or 2. 
Thus, the density is 

ftrfJW, = ^ - 2 "> 2 ' 1 + 6A ' C o" a "" 2 (28) 

v ^ 2n + 3 7 2n + 3 7 - 2 v ; 

If we substitute in the Friedmann equation containing p A , 



a 2 87TG, . Kc 2 

(p + Pa)-—T (29) 



a 2 3 ' a 2 



then we have 



^ = Acga 2 "(g + 3 7 ) i^a 2 - 2 (2-3 7 ) 

a 2 " 3( 2n + 3 7 ) 2n + 3 7 - 2 1 j 

We see that the A term on the right-hand side of eq. (|3~0| ) will fall off 
faster than the curvature and matter density terms so long as c(t) falls off 
fast enough, with 

« < -y (31) 

As expected, this is requires a more rapid increase in c(t) than is required 
to solve the flatness or horizon problems. In particular, in the radiation and 
dust cases it requires c to fall as the Universe expands, with n < —2 and 
n < — 3/2, respectively. 



5 Better Formulations 

The most challenging problem besetting any attempt to examine the cosmo- 
logical consequences of varying c is the formulation of a self-consistent theory 



10 



which incorporates such a variation. The AM approach is not the only way 
of proceeding and so we should look more critically at its relation to other 
things that we know. The most interesting feature in this respect is the fact 
that it permits a description of varying G when c is kept constant. However, 
it is clear that in this case the resulting cosmological equations (f|), © and 
(|]) are not the equations describing the evolution of the Friedmann universe 
in a scalar-tensor gravity theory, like JBD. The reasons are clear. JBD theory 
meets the covariant divergence constraints imposed by the Bianchi identities 
by including the energy density contributed by the scalar field that acts as 
the source of the variations in G. These contributions are not included in the 
AM formulation. In order to improve upon it we could take a scalar-tensor 
gravity theory and carry out the same procedure used to derive (|B|) and check 
that the time variation of c is still permitted. This means that we assume 
that the JBD equations hold for the Friedmann model, with c allowed to be 
a time variable, and then derive the conservation equation from them. For 



JBD theory we have ||, [ I i 



a 

~ 2 

3d 
a 



inp 



!>a to(f) 



Kc 2 (t) 



30 (f>a ' 60 2 
8^ 



(32) 



(3 + 2u)<p 



2 + u)p + 



3{l + u)p 
c\t) 



• 2 

6 low 

7 - -3T> (33) 



+ 



3(26 



in 



3p , 8vrp(4-37) 



, - ■ (34) 

a 3 + 2u" c 2 (t) J 3 + 2lu v ; 

From these three equations we obtain the generalised conservation equation. 
Since after we impose the equation of state, the varying c term only appears 
in the first of these equations, the only new contribution to 2ad is from the 
—2Kcc term. Hence, we have [|TU] 



p + 3- p+4 

a cr 



3Kcd(p 



(35) 



Aira? 

It is instructive to compare this to the AM equation (||]). We see that the 
incorporation of G variation differs. The pG/G term of (|35"D is absent but 
the K term has the same form with G{t) replaced by -1 (£) as usual in JBD 
theory. However, the new system is more constrained because solutions to 
eqns. (|34] ) and ( ^5|) must be consistent. 
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When the matter content is radiation the conservation equation (|33j) has 
the particular solution 

d) = G~ l = constant 



and then fl35|) is identical to the AM equation (^) in the case where G is 
constant but c varies with time. The conditions for resolving the flatness 
and cosmological constant problems then become the same as given for the 
AM formulation above. For other equations of state the situation is more 
complicated. Rearranging the equations, we need to solve 

. 3Kcc<f)a 3 ^ 2 
{pa } = 4~n 

■ = ^ 3 (4-37) 
w } 3 + 2lu 

together with eq. fl32]) . When K = they are solved by the usual exact 
solutions of JBD theory. When the universe is radiation dominated these 
equations have particular solutions with constant which are identical to 
those derived in the last section for the case of varying c and constant G. 



These equations will be studied in more detail elsewhere fll4 |. 

The JBD theory is not the most general scalar-tensor gravity theory. It 
is a particular case of a general class in which the constant BD parameter, 
u, is a function of the scalar field 0, [pTlJl . For these theories, the Friedmann 



universe (with c constant) is governed by the equations (in the A = case) 



■2 



a\ f a\ acj) Kc 2 w(4>)4> 87rp 



2 - + - + 2-^ + 



a J \aj acj) a? 20 0c 2 



a 2 87ip d0 ^(0)0 kc 2 



a 2 30 a0 + 60 2 a 2 ' (36) 

i + 3Hi+ ^ ^pl. ,37, 
2u; + 3 2u + 3 

From these equations we obtain the same generalisation of the mass- 
energy conservation as in JBD (which is the case o;(0) = constant), given by 
eq. (|35|) . Again, radiation-dominated solutions (with = G' 1 constant) and 
varying c of the form found in the last section are always exact particular 
solutions of these equations. 
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6 Some Problems 



6.1 The Velocity Problem 

If c varies here may be a problem with the perturbations to the isotropic 
expansion of the universe which are powers ofv/c. For example, for rotational 
velocity perturbations in matter obeying p = (7 — l)pc 2 , the conservation of 
angular momentum gives pa 4 v = constant. Although the original formulation 
described above does not allow us to impose energy conservation in the usual 
way (E oc c 2 ), we shall assume that angular momentum conservation is 
preserved since it does not incorporate direct c dependence. Hence, we have 

v oc a^-\ 

Now, if 

c(t) = c a n 

we have 

- oc a^~ A ~ n 
c 

If the flatness problem is to be solved, we need 2 — 2n > 37 so we will have 
37 — 4 — n < —2 — 3n. For the A problem to be solved we need 37 + 2n < 0. 
Hence, in the radiation era we need n < — 1 to solve the flatness problem 
and n < — 2 to solve the A problem. Therefore, when 7 = 4/3, we have 

c 

we see that v/c grows in time in both cases as a — > 00. During the dust era, 
7 = 1, and so 

V - oc a— 

c 

where n < —1/2 solves flatness (allowing v/c to decay in —1 < n < —1/2) 
and n < —3/2 solves A (which would always require v/c to grow). 

In general, the condition for solution of the flatness problem has v/c 
growing as a power of a equal to 37 — 4 — n > 9^/2 — 5, while a solution of 
the A problem has it growing as 37 — 4 — n > 9^/2 — 4. So we can only avoid 
blow-up of the velocity perturbations at large a if 7 < 10/9 (if we also want 
to solve flatness) or 7 < 8/9 (if we also want to solve A). 
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6.2 Shear Evolution 



We can extend the AM model to include the effects of simple anisotropics 
with isotropic curvatures. This would correspond to generalisations of the 
Bianchi type I and type V universes in the flat and open cases, respectively 
In the spirit of equations (@)-© the required generalisation which includes 
the shear anisotropy scalar, a, is 

d 2 8nG(t)p a 2 Kc 2 (t) 

a^-^^ + y-^ 2- (38) 

a 4irG(t), 3p 2a 2 

_ = __ii( p+ -^) a __ (39) 

where a(t) is now the geometric-mean expansion scale factor. From these 
equations we may derive the analogue of the matter conservation equation: 

a, p. a , , h . G 3Kcc 

p + 3- (p + ^ + — - (a + 3-a) = -p- + — — . 40 
a & AttG a G 4irGa z 

This equation resembles the conservation equation for two non-interacting 

perfect fluids, one of which (the 'anisotropy' energy density) has the equation 

of state p = pc 2 . However, we see that the original AM prescription leaves 

the material density coupled to the shear evolution. If we add the shear 

evolution equation from the standard (constant c) case then we have 

a + 3-a = (41) 
a 

and the density evolution obeys the same equation as that for the isotropic 
model (a — 0). 

6.3 Black Holes 

It is interesting to note that if a black hole forms with radius 

2GM 

R 9 = — ( 42 ) 

then, from (12) 



Rl = G_ 2 c = _Pa 
R g G c p A 
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This seems to imply that if G is constant and c is falling then primordial 
black holes would grow in size at the same rate as the universe is expanding. 
There will also be consequences for black hole evaporation. The Planck scale 
will grow rapidly oc c~ 5 . It is possible to examine the changes to the laws of 
black hole mechanics using the results for the change in the characteristics of 
black body radiation thermodynamics but we do not digress to discuss this 
here. 



6.4 Late-time fields 

Another aspect of this model is worth commenting on. We have determined 
the rate of variation of c(t) that would allow the flatness (or A) to be solved. 
However, we could also use these equations in the cases where the flatness 
problem was not solved in the same way as a description of a cosmological 
matter field that plays the role of a matter field which can dominate the 
expansion dynamics of the universe at late times. Of course, as with the 
possibility of present-day contributions to the dynamics by the curvature or 
A terms in the case where c is constant, this requires a special tuning of the 
initial sizes of these terms in the Friedmann equation with respect to the 
density term in order that their effects just start to become significant close 
to the present epoch. 

An interesting point to notice about c-varying models is that if the uni- 
verse were to contain a dynamically significant cosmological constant term 
at late times then it would naturally have negligible curvature term and so 
resemble a zero-curvature universe with subcritical total density, just like 
inflationary universes with a residual lambda term. This property allows 



natural almost flat asymptotes to arise, as pointed out in ref. [14] 



6.5 Fine structure constant variations 

One interesting feature of any self-consistent theory of c variation is that 
it enables us to evaluate the consequences of time (and space) variations 
in dimensionless coupling constants like the fine structure constant. In an 
expanding universe with c varying we require pc 2 oc (k B T) 4 /(hcf. Assum- 
ing particle wavelengths are unchanged in the absence of expansion and the 
masses of quantum particles are conserved so h/c remains constant. As the 
universe expands the Planck spectrum remain Planckian but the temperature 
varies as |l[ 
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Hence, 



r -> / i - -2^ 



r 2 

c 2n-l 

a 



for the power-law variation, (|12|)> assumed above. Strong observational limits 
could be placed on any deviation of n from zero using data which establish 
bounds on the temperature of the microwave background at low redshifts. 

The required variation olh/c means that (for constant e 2 ) the fine struc- 
ture constant varies as 

_ e 2 1_ 

he c 2 

If we assume a more modest variation in the evolution of c(t) oc (hit) -1 
then it becomes possible to examine some possible late-time consequences 
of varying a using this model. A detailed particular theoretical formulation 



has been provided by Bekenstein, ||. Elsewhere, Barrow and Magueijo f!3 
have provided a lagrangian-based theory with varying e which differs from 
Bekenstein's (unlike Bekenstein's theory it does not postulate that the gravi- 
tational field equations remain unchanged). This theory may be transformed 
by a suitable change of units into a theory with constant electron charge e 
and varying c like the one described above, or it may be presented as a theory 
with constant c and varying e. This features displays explicitly that invariant 
meaning can only be ascribed to variations of dimensionless 'constants' (in 
this case the fine structure constant). 

A more far-reaching aspect of c-variation is the fact that a fall in the 
value of c at very early times would lead to a strengthening of dimensionless 
gauge couplings unless other changes in the structure of physics could offset 
this. Thus, attempts to reconstruct the history of the very early universe 
could not count on the applicability of asymptotic freedom and the ideal gas 
condition at early times in the usual way. 

7 Discussion 

The contemplation of a variation in the speed of light in the very early 
universe presents a host of self-consistency problems and boundary effects 
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with other parts of the core of modern physical theory. Unlike the case of 
varying G, it requires deep structural changes to many of the foundations 
of physics. In this paper we have considered a minimalist varying-c theory 
recently proposed by Albrecht and Magueijo [|] as a new way of solving 
the horizon and flatness problems. However, unlike inflation, it also offers 
a way of solving the cosmological constant problem. Whereas Albrecht and 
Magueijo modelled changes in c as a sudden fall in the universal value of 
c,occurring as if at a phase transition, we have considered the behaviour of 
variations which vary as a power of the cosmological scale factor in order to 
determine the rates of change that are needed to solve different cosmological 
problems. Our main results are as follows: 

a. If the velocity of light varies with the cosmological scale factor, a(t), 
as c = coa n and the equation of state of matter is p = (7 — l)pc 2 , then 
the curvature term becomes negligible for the expansion of the Universe as 
a — > 00 if c falls fast enough, with n < ^(2 — 37), regardless of the behaviour 
of G(t). The same condition allows for a resolution of the horizon problem. 
This generalises the conditions for the inflationary resolution of the flatness 
and horizon problems which apply when c is constant (n = 0) and G is 
constant. 

b. If the velocity of light varies with the cosmological scale factor, a(t), as 
c = Coa n ,the gravitation 'constant' varies as G(t) = Goa g , and the equation 
of state of matter is p = (7 — l)pc 2 , then the condition for a vacuum stress 
with pa — ~Pa° 2 i associated with a geometrical cosmological constant term 
to become negligible as a —>■ 00 is that c fall at a rate such that n < —3^/2. 

c. We have shown how it possible to improve the formulation of the 
problem in order to recover a scalar-tensor theory describing G variation 
in the limit of constant c. This formulation only allows the solutions of 
the AM theory to persist in the case when the equation of state is that of 
black body radiation. This problem also suggests that we should attempt a 
fuller formulation of the varying-c theory by deriving the variation of c from 
some scalar field that contributes an energy density to the evolution of the 
Universe. 

d. We have highlighted some potential problems with the evolution of 
rotational velocities as the universe expands. It is also worth noting that 
the varying-c model does not solve the flatness and horizon problems by 
means of a period of evolution that is close to that described by a time- 
translation invariant de Sitter space-time. As a result there does not appear 
to be any distinctive set of fluctuations that should emerge from a period of 
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cosmological evolution in which c changes. 

e. Possible variations in c lead to changes in the fine structure constant 
and to other gauge couplings in the very early universe if no other changes 
in physics exist. However, if the variation in c is confined to a very early 
period of evolution soon after the expansion commenced there need be no 
conflict with astronomical constraints on any time or space variation in the 
fine structure constant. 

f. A varying-c theory can be transformed in to a varying-e theory with 
constant c by a suitable transformation of units. Theories of this sort can be 
given a lagrangian formulation if required. 

We believe that these features of the naive varying-c model proposed by 
Albrecht and Magueijo are sufficiently interesting for it to be worthwhile 
exploring a fuller, more rigorous formulation of a varying-c gravity theory 
based upon an action principle which would enable its consequences for other 
aspects of physics and cosmology to investigated more rigorously. 
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